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Abstract 

Recent claims of the Seattle group of evidence of "slope difference between CD [Coulomb Dis- 
sociation] and direct [capture] results" are based on wrong and selective data. When the RIKEN2 
data are included correctly, and previously published Direct Capture (DC) data are also included, 
we observe only a 1.9 sigma difference in the extracted so called "scale independent slope (b)", 
considerably smaller than claimed by the Seattle group. The very parameterization used by the 
Seattle group to extract the so called b-slope parameter has no physical foundation. Considering 
the physical slope (S' = dS/dE) above 300 keV, we observe a 1.0 sigma agreement between slopes 
(S') measured in CD and DC, refuting the need for new theoretical investigation. The claim that 
Sn(0) values extracted from CD data are approximately 10% lower than DC results, is based on 
misunderstanding of the CD method. Considering all of the published CD £17(0) results, with 
adding back an unconfirmed E2 correction of the MSU data, yields very consistent £'17(0) results 
that agree with recent DC measurements of the Seattle and Weizmann groups. The recent cor- 
rection of the b-slope parameter (0.25 MeV -1 ) suggested by Esbensen, Bertsch and Snover was 
applied to the wrong b-slope parameter calculated by the Seattle group. When considering the cor- 
rect slope of the RIKEN2 data, this correction in fact leads to a very small b-slope parameter (0.14 
MeV" 1 ), less than half the central value observed for DC data, refuting the need to correct the 
RIKEN2 data. In particular it confirms that the E2 contribution in the RIKEN2 data is negligible. 
The dispersion of measured £17(0) is mostly due to disagreement among individual DC experiments 
and not due to either experimental or theoretical aspects of CD. Additional uncertainty exists due 
to theoretical extrapolation procedure that for example was recently estimated by Descouvemont 
to be 6%, more than a factor of two larger than suggested by the Seattle group. Uncertainty of 
the slope (S'), that thus far was not measured with high precision, leads to a substantial error of 
extrapolated 617(0). 

PACS numbers: 26.30.+k, 21.10.-k, 26.50.+k, 25.40.Lw 



*Work Supported by USDOE Grant No. DE-FG02-94ER40870. 
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I. INTRODUCTION 



The method of Coulomb Dissociation (CD) was developed in the pioneeringwork of Baur, 
Bertulani and Rebel |l| and has been applied to the case of the CD of 8 B . These 

data were analyzed with a remarkable success using only first order Coulomb interaction 
that includes only El contribution. Indeed early attempts to refute this analysis by intro- 
ducing a non-negligible E2 contribution were shown J(| to arise from a misunderstanding of 
the RIKEN1 data; hereafter all collaborative experiments are identified by the location of 
the laboratory where they were performed. Later claims for evidence by the MSU group []| 
of non-negligible E2 contribution in inclusive measurement of an asymmetry, were not con- 
firmed in a recent exclusive measurement of a similar asymmetry by the GSI2 collaboration 

It was then a great surprise to learn in a recent Physical Review Letter by Esbensen, 
Bertsch and Snover 8] that higher order terms and an E2 contribution are an important 



correction to the RIKEN2 data |8j . They quote in their PRL paper [8] statements from the 
Seattle group l| and indeed the PRL paper is based on these claims We demonstrate 
that these claims are based on misrepresentation of data and a misunderstanding of the CD 
method. 



II. MISREPRESENTATION OF DATA ON B-SLOPE PARAMETER 

The Seattle group chose to parameterize the astrophysical cross section factor, Sn, as S = 
a(l+bE), and they show in Fig. 19 of their paper [9j the so called "scale independent slopes 
(b)" extracted for energies smaller than 425 keV and between 830 - 1300 keV. We first note 
that excluding the region of the 632 keV resonance is possibly relevant (but not necessary) 
for Direct Capture (DC) data, but in CD data the contribution of that state is reduced 
(by a factor of j3 2 ) and is negligible. For the RIKEN kinematics it is in fact smaller than 
the quoted error bar. The GSI CD data are quoted with that small contribution already 
removed and there is no reason to exclude a certain energy region to determine the slope of 
CD data. 

In Fig. 1 we show the b-slope parameter extracted from both DC and CD data in the 
energy range of 300 - 1400 keV. For the DC data we exclude the same energy region <| 
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FIG. 1: The extracted so called "scale independent slopes (b)" of world data. We correctly 



plot 



the RIKEN2 data and add several data points that were neglected by the Seattle group £ 
discussed in the text. The range of "average values" is indicated and discussed in the text. 



as 



due to 632 keV resonance and we subtract its Ml contributions at higher energies. All so 
obtained fits have an acceptable x 2 / u close to unity. The obtained b-slope parameters agree 
with Ref. 0, however the so called scale independent b-slope parameter of the RIKEN2 
data 3] was plotted incorrectly by the Seattle group in Fig. 19 of their paper It is 
plotted correctly in Fig. 1, where we also include three DC data points that were omitted 
by the Seattle group |9j, as we discuss below. 

A detailed discussion of the difference between the RIKEN2 and GSI1 data was already 
presented in the literature fl. For exatnple f r0m F ig . 4 of the GSI2 paper Q it i 8 dear that 
the central value of the slope of the RIKEN2 data is smaller than the central value of the 
GSO data D but they ate plotted by the Seattle group in Pi, 19 with almost the sa,„e 
central value. We can not explain the origin of this simple mistake of the Seattle group 
[l^ . which in of itself might be of no consequence, except for its relevance for the newly 
published PRL paper Nl as we discuss below. 

nn 

Very recently it was also suggested [lfj, |11[ that a reanalysis of the GSI2 data may yield 
a smaller slope for the GSI2 data, as also shown in Fig. 1. This latest analysis has yet to be 
scrutinized p| and published, none-the-less it is shown here since it was already presented 



in international meetings 
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□ 



111 ] . But this result |l0j, [ll( , as well as all other non published 



results, are not included in the calculation of the average slope and extrapolated Si7(0) 
values of CD data. 

The published RIKEN2 Q, GSIlQ, GSI2 Q, and MSU data Q, yield a 1/a weighted 
average for the "b-slope parameter", b = 0.51 ±0.06 MeV~ x with x 2 jv= 1.1. With a central 
value that is slightly smaller than quoted in Ref. |9j, b = 0.55 ±0.08 and x 2 /V= 0.2. The 
quoted small x 2 / u is due to the wrong value used for the RIKEN2 data in Ref. We 
also note with much disappointment that the discussed average b-slope parameters jfl| are a 
factor of 10 larger than the very same b-slope parameters shown in Fig. 19 Q]. The smaller 
values plotted in Fig. 19 are the correct ones. 

In the same figure the Seattle group suggests that a remarkable agreement exists between 
_d b- S ,ope parsers of DC data Q Q Q Q. The DC data shown in Efc 
19 [9J] does not include all available data. We note that the "old" data of Filippone (16J 
are included [^], but the data measured by Vaughn et al. Parker [3], and Kavangh 

et al. 3] were ignored in their Fig. 19. The last three data sets were deemed j^] not 
useful at low energies for extrapolating Si7(0), but they are certainly useful for studying the 
slope of the data measured at energies between 300 and 1400 keV as done in Fig. 19 of the 
Seattle group. Most disturbing is the fact that they exclude the published data of Parker 
Q| and Vaughn et al. but included that of Filippone ^(|. As we show below including 
these data as well as including the correct value of the RIKEN2 data reduces the claimed 
disagreement by a factor 2. 



III. DISAGREEMENT AMONG DIRECT CAPTURE DATA 



It is well known that the "old" data on DC 7 Be(p,'j) 8 B reaction [ifij, lift |l2j] exhibit 
major systematic disagreements. But the situation is not improved with "modern" data on 
DC 0,13, 14, 15 1 ■ The data of the Orsay group [rJ and Bochum group [l^ do not agree 
with that of the Seattle group and Weizmann group [l5 1 . The disagreement is by as much 
as five sigma and there is not a single measured data point of the Bochum group that agrees 
with a data point measured by the Seattle group. Most disturbing is the disagreement of 
the b-slope parameter extracted from the two "modern" high precision measurements of the 
Seattle and Weizmann jl^ groups, shown in Fig. 19 of as well as in Fig. 1 of this 
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japer. From Fig. 19 of |l| it is clear that the b-slope parameter of the Orsay 13], Seattle 
J and Weizmann data do not have overlapping error bars. As we show below the, 
disagreement between the slope parameters is in fact more significant. The slightly better 
agreement between "modern" DC data is an artifact of the parametrization used by the 
Seattle group. 

The large systematical disagreement between DC data can not be handled algebraically 
to extract a meamngful average slope For examp.e the slopes of the Orsay Q, Seattle fl 
and Weizmann [15j shown in Fig. 19 |9J can not be all true if the error bars shown in Fig. 19 
[^| are correct. None-the-less it has been customary to artificially enlarge the error bars by 
multiplying it by the square root of x 2 1 v i so as to make data with systematical differences 
appear as if it is statistically distributed. Using such a procedure for all published DC 
data Q M, 15, 16, |l7, |l8| and the BE1 and BE3 data sets of the Seattle group Q 



we 



extract a 1/cr weighted average b-slope parameter, b = 0.38 ±0.02 MeV~ l with x 2 l v= 1-8- 
Taking into account the bad x 2 l v as discussed above, we obtain b = 0.38 ±0.03 MeV~ x as 
compared to b = 0.311 ±0.014 MeV~ l and x 2 l v= 1-9, quoted by the Seattle group. 

We conclude that the b-slope parameter can not be extracted from DC data with the (im- 
pressive) accuracy of 4.5% |9j, unless one excludes some of the DC measurements discussed 
above. An error which is approximately a factor of 2 larger seems like a more reasonable 
choice. Also a difference of the extracted b-slope parameter of 0.25 MeV" 1 quoted by 
the Seattle group is not confirmed in this analysis where we find a difference which is 
approximately a factor 2 smaller. 

The so obtained average b-slope parameters in CD, b = 0.51 ±0.06 MeV -1 , differs 
by only 1.9 sigma from that extracted for DC data, b = 0.38 ±0.03 MeV -1 , and can 
hardly be characterized as significant. This slight difference is particularly negligible in 
view of the substantial systematical disagreement among individual data sets measured 
in DC experiments. Furthermore, it should not be considered sufficient to motivate new 
theoretical investigation or for that matter a publication in the prestigious Journal of the 
Physical Review Letters. As we show below the actual differences of measured slopes are 
even smaller. 
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IV. LACK OF PHYSICAL JUSTIFICATION OF THE SEATTLE SLOPE PARAM- 
ETERIZATION 



We now demonstrate that the parametrization S = a(l + bE) used by the Seattle group 



01 



has no physical justification. We refer the reader to Baye's seminal paper |2l( on the theory 
of DC where we find that for an external capture reaction (and only when the conditions 
for an external capture reaction are satisfied, i.e. below 100 keV): 

S(E) = 5(0) [1 + Sl xE] (equ 1) 
5(0) = 5,(0) + 5,(0) (equ 2) 



and 



si = §f[si s + 8 U x |M] (equ 3) 



Equation 1-3 were derived for external capture and are correct only when the conditions 
for external capture are met. Specifically the value of S\ is explicitly negative as predicted 
for energies below 100 keV. In such a case the logarithmic derivative S'/S(0) is shown to be 
an invariant j^j. 

At energies above 300 keV the observed slope are manifestly positive. For higher energies 
one must add higher orders beyond the linear term used in equ 1, so as to turn the sign 
of the slope from negative to positive. For example the Taylor expansion of the theoretical 
curve predicted by Descouvemont and Baye 22] can not be truncated below third order 
to yield a reasonable representation of the predicted curve at all energies up to 1.5 MeV. 
Clearly at energies above 300 keV the truncation of the Taylor expansion to a linear term, 
leads to unphysical expansion. 



Furthermore, as we show in 
constant as a function of energy 



? ig . 2 the slope of the d-wave component is essentially 
231 ] (approximately +10 eV-b/MeV), but the slope of the 
s-wave component varies with energy between approximately -20 eV-b/MeV at zero energy 
to approximately -3 eV-b/MeV at higher energies. The overall slope (S' = dS/dE) is a sum 
of the two components which is energy dependent. At very low energies the s-wave negative 
slope is dominant and is almost constant hence the linear truncation as in equ. 1. At 
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energies above 500 keV the d-wave slope dominates and the variation of the s-wave slope is 
small leading to almost constant positive slope. At high energies the linear dependent of the 
S-factor with a positive slope is an artifact of the d/s ratio. But the over all slope at higher 
energies is dependent on the low energy cut of the data. In particular the fit parameters 
"a" used in |9| is most sensitive to the choice of the low energy cutoff of the data between 
300 and 500 keV, and the so-called scale independent b-slope parameter varies due to the 
selected range of data and not due to the physical slope (S'). 

In equ. 1 the over all normalization factor S(0) is directly related to the astrophysical 
cross section factor at zero energy, as well as the ANC of the physical wave function j^j]. 
But the fit parameter "a" used by the Seattle group [l| has no physical meaning and is not 
directly related to the ANC or the spectroscopic factor. 

Furthermore, the observed positive slope (S' = dS/dE) of the cross section factor at 



energies above 500 keV is direc 



try re 



Q, 21, Q, 24, Q, Q 



ated to the (model dependent) d-wave contribution 



2J, |25j, 26]. It is self evident that the d/s ratio is 



designated by d/s ratio 
independent of S(0), since both d and s components are proportional to the same ANC [2lj. 
Hence it is clear that at energies above 300 keV the slope (S') is not directly related to the 
over all normalization S(0) and the so called "scale independent slope (b)" defined by the 
Seattle group 9] has no physical justification. 

To illustrate thispoint we show in Fig. 2 the s and d waves contributions calculated 
by Jennings et al. [23]. In the same figure we show the predicted S-factor with 5*^(0) that 
is increased from 1.3 eV-b 23] to 1.7 eV-b, which is certainly within the limit of accuracy 
of theoretical predictions. This yield a very insignificant (1.8%) change of S(0) but a very 
significant change (~ 25%) of the observed slope (S') above 300 keV. This schematic model 
most vividly demonstrates the fact that the parametrization S = a(l+bE) has no physical 
justification nor does the assumption that S'/S(0) above 300 keV is an invariant. 

We add that a knowledge of the d-wave component at measured energies (e.g. above 300 
keV) is essential for extrapolating 5*17(0). For example at 500 keV the d-wave contribution 
amounts to 30% of the measured Sn and above 1,150 keV it is dominant. At zero energy, 
on the other hand, the d-wave contribution is predicted to be small (Sd(0)/S s (0) ~ 6%). 
Thus in order to accurately extrapolate Sn to zero energies one must remove the d-wave 



contribution from measured 5*17, as has been emphasized long ago by Robertson 



2J. This so 



far has been done by means of theoretical estimate of the d-wave component and a chi-square 
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FIG. 2: The s and d wave contributions calculated by Jennings et al. [23( and with an increased 
d/s(0) = 8.5% value. 

fit of data by the predicted s + d wave components. 

But currently there is no direct way to test the validity of the model dependent prediction 
of the d-wave component. In fact the slope (S') seems thus far the only way to determine 
the d-wave contribution. As we show below the slope has not been measured with high 
accuracy and large discrepancies still exist between measured DC data. Indeed the ill defined 
theoretical d/s ratio leads to an uncertainty of the extrapolation. Other theoretical issues 
were also discussed by Descouvemnot and they yield an uncertainty of at least 6% j^, a 
factor of 2 larger than the theoretical uncertainty quoted by in Ref. 

V. THE SLOPE S' 

For that reason we plot in Fig. 3 the extracted slope S' = dS/dE for data measured 
between 300 and 1400 keV, with the same exclusion of DC data due to 632 keV resonance 
between 450 and 850 keV, and correcting for the Ml contribution of that resonance above 
850 KeV as discussed in section II. Figure 3 includes all published data and the slope 
of the RIKEN2 data is plotted correctly. In this figure we note a sharper disagreement 
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FIG. 3: The measured slopes (S' = dS/dE) of world data measured between 300 and 1400 keV, 
as discussed in the text. The range of "average values" is indicated and discussed in the text. 

between "modern" DC data. In particular we extract three different slopes values for the 



13. 



141 ]. an average of the 



four "modern" experiments on DC: a Bochum-Orsay low value 

three Seattle results |9] and the Weizmann higher value. [3]. These three values of "modern" 
DC slopes do not exhibit overlapping error bars. In sharp contrast all four published CD 
slopes shown in Fig. 3 are in agreement within the quoted error bars. 

From these data we extract for CD data the l/a weighted average S' = 7.9 ±1.0 eV- 
b/MeV with x 2 /v= 0.5 and for DC data S' = 6.7 ±0.5 eV-b/MeV with x 2 /"= 1-5, hence 
S' = 6.7 ±0.6 eV-b/MeV. Clearly the average slope extracted in CD data agrees (within 
1.0 sigma) with that extracted from DC data. This agreement is considerably better than 
observed between individual DC measurements, as we discuss above. 

We conclude that there is no significant disagreement between slopes extracted from CD 
and DC data. This removes the very raison d'etre for the recent PRL paper 



10 



VI. COMPARISON OF RIKEN2 DATA WITH THEORY OF ESBENSEN, 
BERTSCH AND SNOVER 

Furthermore, in this PRL paper [8] we find a substantial (50%) correction of the b-slope 
parameter extracted for the RIKEN2 data that is implied to be 0.25 MeV" 1 . As we discuss 
above the b-slope parameter was computed incorrectly by the Seattle group, as shown in 
their Fig. 19 |<g. For the RIKEN2 data b = 0.39 ±0.08 MeV' 1 , and the said correction 
yields: b = 0.14 ±0.08 MeV" 1 which is considerably (more than a factor of 2) smaller 
than the observed central value for DC data, b = 0.38 ±0.03 MeV~ l . We conclude that 
;he corrections suggested by Esbensen, Bertsch and Snover, in sharp contrast to the claim 
^], in fact lead to a disagreement between the RIKEN2 data and DC data, and are not 
relevant at least for the RIKEN2 data, and in particular they once again confirm that the 
E2 correction is negligible. 



VII. Sn(fi) EXTRACTED FROM CD DATA 
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In Fig. 20 of the Seattle paper [9j] they show extracted 617(0) from CD using the extrap 



olation procedure of Descouvemont and Baye |22j, and based on this analysis it is stated 
[^] that "the zero-energy extrapolated 617(0) values inferred from CD measurements are, 
on the average 10% lower than the mean of modern direct measurements". The extracted 
617(0) shown in Fig. 20 [9j are only from data measured at energies below 425 keV and the 
majority of CD data points that were measured above 425 keV were excluded in Fig. 20 
This arbitrary exclusion of (CD) data above 425 keV has no physical justification. For 



example as shown by Descouvemont 



271 ] the theoretical error increases to approximately 



5% at 500 keV and in fact it is slightly decreased up to approximately 1.0 MeV, and there 
is no theoretical justification for excluding data between 500 keV and 1.0 MeV. In DC 
measurements the well known contribution of the 632 keV resonance needs to be subtracted, 
but that is not the case for example in RIKEN CD data. Furthermore, as we discussed above, 
the slope of the data between 300 and 1,400 keV is essential for determining the d/s ratio 
and the extrapolation to zero energy. Excluding data above 425 keV reduces our sensitivity 
for testing the various model prediction of the d/s ratio. 

When including CD data measured only at energies below 425 one runs into a more serious 
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systematical problem. Namely, the relative energy measured in CD (E re i) is determined 
mostly from the proton- 7 5e relative angle (#17). At small relative energies (as well as at 
small scattering angles, 9g) plural scattering in the target are of major concern. These are 
estimated theoretically and are known to be inaccurate. The effect of plural scattering in 
the target is indeed known for practitioners in the field of CD and has been discussed on 
several occasions. It leads to a systematical uncertainty of the measured Sn of the order of 
2 eV-b. 

Furthermore, the yield of the CD of 8 B arises from a convolution of the nuclear cross 
section which is rapidly dropping toward low energies, and the virtual photon flux that is 
rapidly increasing toward low energies. This generate a yield that is almost constant (±20%) 
at E re i = 300 - 800 keV. Note that over the same energy range the DC yield changes by 
almost a factor of 10. 

Thus when excluding the CD data above 425 keV, the Seattle group excluded the data 
that were measured with the best accuracy and with smallest systematical uncertainty. This 
exclusion of "the best energy region" measured in CD experiment is a manifestation of their 
misunderstanding of the CD method. If in fact one insists on such a flawed analysis of CD 
data, one must estimate the systematic uncertainty due to this selection of data. This has 
not been done in the Seattle re-analyses of CD data P|. 

Instead we rely here on the original analyses of the authors. In Fig. 4 we show the 
617(0) factors extracted by the original authors who performed the CD experiments. These 
results include all measured data points, and are consistently analyzed with the extrapolation 
procedure of Descouvemont and Baye J3] . The potential model of Typel [s| was also used 
in the GSI2 paper, but the so quoted Si7(0) are not shown in Fig. 4. 

We note that the five experiments on the CD of S B show a remarkably good 

agreement within the quoted error bars, in sharp contrast to the confusion that exists in 
"old" 0, 0,0,0 and "modern" DC results 0,0, QQ. The results of the RIKEN- 
GSI experiments must be considered as a continuation of the same experiment (essentially 
the same experiment repeated four times) with an improved kinematical and experimental 
conditions. Thus the four results can not all be assigned the same weight. 

We note that while these four results are consistent within the quoted error bar, they 
show a systematic trend of an increased Si7(0) (to approximately 20.7 eV-b), while the 
error bars are reduced. The MSU result on the other hand includes a model dependent E2 
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FIG. 4: Measured Si7(0) as originally published by the authors who performed the CD exper- 
iments. These analyses include all measured data points , lZ| using the extrapolation 
procedure of Descouvemont and Baye [24]. We also plot the MSU data with the E2 correction 
(~ 8%) [?J added back to the quoted 5*17(0), as discussed in the text. The range of Si7(0) results 
from the measurements of DC by the Seattle jj| and Weizmann groups is indicated. 



which was not confirmed in a 



correction 8%) deduced from inclusive experiments 

recent exclusive measurement of a similar asymmetry When this E2 correction is added 
)ack to the quoted MSU result 7], as shown in Fig. 4, together with the published RIKEN2 
jj, GSI1 Q, and GSI2 Q results we obtain a 1/a weighted average of S 17 (0) = 20.0 ±0.7 

with y 2 / v = 0.5, which is in excellent agreement with the measurement of the Weizmann 



group 



15| and in agreement with the measurement of the Seattle group 



□ 



We conclude that the statements (which were repeated in the literature) on substantial 
differences between Si 7 data from DC and CD measurements as well as extracted £17(0) are 
based on misrepresentation of data, misunderstanding of theory and a biassed approach of 
the Seattle group Quite to the contrary we find a good agreement among CD measure- 
ments as well greement of the CD results with the two most recent high precision 
results of DC experiments. Further attention must be given to an accurate measurement of 
the slope (S') and the d-wave contribution to 6*17 measured at energies above 300 keV so as 
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to allow accurate extrapolation to zero energy. 
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